Analysis of Low Energy Pion Spectra 
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The transverse mass spectra and the rapidity distributions of 7r + and tv~ in Au-Au collisions at 
2, 4, 6, and 8 GeV-A by E895 collaboration are fitted using an elliptically expanding fireball model 
with the contribution from the resonance decays and the final state Coulomb interaction. The 
ratio of the total number of produced tt~ and tv + is used to fit the data. The resulting freeze-out 
temperature is rather low(T/ < 60 MeV) with large transverse flow and thus resonance contribution 
is very small. The difference in the shape of mt spectra of the oppositely charged pions are found 
to be due to the Coulomb interaction of the pions. 

PACS numbers: 24.10.Pa,25.75.-q 



Pion production just above the threshold energy is 
quite different from that at very high energies such as 
RHIC energy since the ratio of ir~ to 7r + at very high en- 
ergies is one which is not the case at low energies. At just 
above the threshold energy, pions are produced through 
the production of A resonances and counting all the pos- 
sible channels of A decay the difference in the compo- 
sition of isospins in the colliding nuclei appears as the 
difference in the numbers of the two oppositely charged 
pions [J, 0, El, whereas at high energies many channels 
producing pions are open and small asymmetry in the 
initial isospin does not matter. 

Other features of the pion spectra at low energies 
are [J H 0, 13, & @, 0]: (l)The transverse momentum 
spectra both of the n~ and tt + seem to have two temper- 
atures. Usually the low temperature component in the 
low momentum region is attributed to the pions decayed 
from resonances, especially the delta resonance, while the 
higher temperature component in the mid-momentum re- 
gion is the thermal ones. (2) Transverse momentum spec- 
tra of ir~ and ir + at very small momentum are different 
in the sense that the 7r + spectra is convex in its shape 
while the ir~ spectra does not show this behavoir. This 
difference in low momentum region is due to the Coulomb 
effect. The hadronic matter formed during the collision 
has charge which comes from the initially colliding two 
nuclei and thus the thermal pions escaping from the sys- 
tem experience the Coulomb interaction. The Coulomb 
interaction of 7r~ may bend the spectrum in the low mo- 
mentum region upward and thus it is hard to disentan- 
gle the contribution from the delta resonance and the 
Coulomb interaction in the low momentum region. (3) 
Width of the rapidity spectra of ir~ and tt + are much 
wider than those from the isotropically expanding ther- 
mal model. The wide width may either come from partly 
transparent nature of the collision dynamics or the ellip- 
soidal expansion geometry. In order to fit large rapidity 
width using expanding fireball model one usually needs 
large longitudinal expansion velocity. 



Even though all those features mentioned above are not 
new, calculations with all those features put in together 
comparing each contributions in detail is hard to find. 
There are claims that the properties of A resonance are 
modified inside the hadronic matter formed even at this 
low energy 5. 6]. In order to draw any conclusion, one 
should have a model which can explain all of the features 
above mentioned. Lacking one or two features in the 
calculations, the result may not be conclusive. 

In this paper, we analyze the pion spectra in Au+Au 
collisions at 2,4, 6, and 8 A-GeV measured by the E895 
collaboration^ using the expanding fireball[l, 0, [13] 



and the final state 
The geometry of 
and can be varied to 
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with the resonance contribution 
Coulomb interaction P, [TJ, 0j|, 
the expansion used is ellipsoidal 
sphere and cylinder, by taking the transverse size R as a 
function of the longitudinal coordinate, z. 

At just above the threshold energy, pions are produced 
through the production of A resonances and their sub- 
sequent decays. The ratio of 7r~ and tt + is given from 
the initial isospin conservation as ^ = g^a+j^f ~ 1.94 
for Au+Au collisions 0, 0, Q . Hence it is expected that 
at low beam energies near 2 GcV-A in Au+Au collisions, 
the ratio is near 1.94 and then as the beam energy is in- 
creased the ratio will decrease eventually to one. In the 
present calculation, the ratio of normalization constants 
for 7r~ and 7r + , R is taken as a fit parameter in order to 
investigate the beam energy dependence of R. 

We assume that once the pions are produced they 
rescatter among themselves and thermalize before they 
decouple from the system. Hence we assume thermaliza- 
tion of the pionic matter. However, the total number of 
negatively and positively charged pions are not in chem- 
ical equilibrium and the ratio is governed by the isospin 
asymmetry of the initially colliding nuclei. We keep the 
ratio as a fitting parameter and want to compare with the 
expected value of 1.94 near threshold of the pion produc- 
tion. 

After the formation of a pionic fireball, it expands and 
cools down until freeze-out when the particle production 
is described from the formalism of Cooper- Frye[l4j. For 
the equilibrium distribution function we use the Lorentz- 
boosted Boltzmann distribution function, where is the 
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TABLE I: Fitted values for each parameters. 



E895, Au+Au at 2A-GeV 
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expansion velocity the space-time of the system. 

E^L = f p »da^x)f(x,p) (1) 



where 



d 3 P (2tt)3 J Ej 



f{x,p) = exp( - ) 



(2) 



For the geometry of the expanding fireball[8|, we as- 
sume azimuthal symmetry and further assume boost- 
invariance collective dynamics along the longitudinal 
direction [l5j. In this case it is convenient to use as the 
coordinates the longitudinal proper time r = y/t 2 — z 2 , 
space-time rapidity rj = tanh _1 (z/i) and the trans- 
verse coordinate r±. Then the 4- velocity of expansion 
can be expressed as u^(x) = 7(1, v±(x)e r , v z (x)) where 
V z (t, rj_,ry) = tanh?7, which is the result of the longitudi- 
nal boost-invariance. In the transverse direction we take 
a linear flow rapidity profile, tanh -1 i^ = p(r))(r±/ Rq) 
where i?o is the transverse radius at midrapidity. Here 
one takes p(rj) — po\A — (f) 2 /lLJ f° r tnc elliptic ge- 
ometry and a constant value of p, i.e. p(r/) = p f° r the 
cylindrical geometry. As is the same for the SPS energy 
by Dobbler et. at, the elliptic case fits the pion spectra 
a little better. 

As pions escape from the system at freeze-out, they ex- 
perience the Coulomb interaction with the charge of the 
system which are mainly due to the initial protons in the 
colliding nuclei. The Coulomb effect on the particle spec- 
tra are studied in detail in refs.[ll|, E2, EH for the static 
and dynamical cases. Here due to the low beam energy 
we restrict ourselves only to the static case. The system 
is expanding rapidly in the longitudinal direction and the 
change in the longitudinal momentum is negligible. Only 
the transverse momentum of the charged particles will be 
shifted by an average amount [11 



p c = Ap A 



2e z 



,dN ch 1 



dy R 



f 



(3) 



where Rf is the transverse radius of the system at freeze- 
out. Due to the lack of detailed knowledge, p c is taken 
as a fit parameter in the present calculation. And in this 
way the beam energy dependence of p c can be studied. 

Since the transverse momentum of the escaping ther- 
mal pions are shifted by the amount p c , i.e. p t = pt t o±p c . 
the invariant cross section can be written in terms of the 
unshifted momentum (p t .o,yo)- 
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FIG. 1: Rapidity distribution of 7r and ir + in Au+Au col- 
lision at 2 A-GeV by E895 collaboration^. 



d 3 N d 3 N dp 3 E 



(4) 



where (E^f)a is the unshifted invariant cross section. 
Integrating over the rapidity y one gets the equation for 
the transverse momentum spectra, 



dN I dy d 2 N Ptfi 
Ptdpt J Pt,odpt,ody p t 



(5) 



and the rapidity spectra is obtained by integrating over 
the transverse mass. 



dN 
dy 



mtfidmtfii 



d 2 N pt |0 
Ptfidptfi p t 



(6) 



Finally one has to add the contribution from resonance 
decay[10| to both the transverse spectrum and rapidity 
distribution. Here we assume that the resonances decay 
far outside the system and the Coulomb interaction of 
the pions decayed from the resonances is neglected. 

The fitted values for the parameters are tabulated in 
the Tab. 1 and the results of the fitting are shown in 
Figs.1-2 for 2 A-GeV. The fitted value for the ratio R = 
7r~ /tt + is close to 1.94 at 2 and 4 A-GeV as expected and 
decreases to 1.38 at 8 A-GeV, which eventually becomes 1 
at higher energies such as RHIC energies. In other words, 
at this very low energy pion isospin is not in chemical 
equilibrium. 
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FIG. 2: 7T~ transverse mass spectrum for each rapidity 
bin(Aj/ = 1.0) measured by E895 collaboration 0] in Au+Au 
collisions at 2 A-GeV. Data on the top line is for the ra- 
pidity bin —0.65 < y < —0.55 and the next one is for 
-0.55 < y < -0.45 scaled by 0.1, etc. 

The freeze-out temperature is rather small, Tf < 60 
MeV and the expansion velocities in both the longitudi- 
nal and transverse direction are quite large. The large 
longitudinal expansion velocity(> 0.8c) is needed to fit 
the large width of the rapidity distribution. The low 
freeze-out temperature together with the large transverse 
expansion fits the transverse spectra of pions quite nicely. 
If not for the large expansion velocity, one usually gets 
much larger freeze-out temperature(Ty > 80 MeV. 

Since the freeze-out temperature is small, there are 
very few resonances at freeze-out, especially A, and thus 
the resonance contribution is negligible. This is reason- 
able since at this low energy near the pion threshold en- 
ergy, production of particles with mass larger than pions 
is rare and their contribution to the pion spectrum is 



negligible. 

The pion transverse momentum spectra looks like that 
there are two slopes; one for small momentum region and 
another for the higher momentum region near the pion 
mass. The smaller slope at lower transverse momentum 
is usually attributed to the pions from the resonance de- 
cay, especially from the A decay [HQ. However, present 
calculation shows that this is not the case at low beam 
energies. 

The shape of the transverse mass spectra of tt~ and 
7r + are different especially in the small mass region. As 
the transverse mass m t decreases, the m t spectrum of n~ 
increases sharply while that of ir + saturates showing the 
convex shape. This difference is due to the Coulomb in- 
teraction of pions leaving the system which has the charge 
from the initially bombarding nucleons. The change in 
the transverse momentum due to the Coulomb interac- 
tion decreases from 25 GeV/c at 2 A-GeV to 15 GeV/c at 
8 A-GeV. This behavior can be understood from the in- 
crease of the screening effect since the number of charged 
pions increase at higher energies. At very high energies 
such as RHIC energies, the momentum change from the 
Coulomb interaction will be small. 

The emerging picture of pion production at low en- 
ergy is that the pions are produced through the inter- 
mediate A formation and thus they are not in chemical 
equilibrium in isospin. They make collisions and ther- 
malize to form a fireball which expands and cools until 
the freeze-out. Since the fireball has charge which is from 
the initially colliding nucleus, the pions leaving the sys- 
tem experiences the Coulomb interaction which makes 
the difference of the m t spectra of the two oppositely 
charged pions. 
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